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In linear response theory, the dielectric 
response at zero frequency sometimes 
appears to violate the /-sum rule, which 
has apparent implications for causality. 
Here, we study the origin of this apparent 
discrepancy, focusing on Lindhard's 
formula for the transverse response of 
the electron gas. At non-zero frequency, 
first-order poles contribute to the 
imaginary part of the dielectric function 
in the usual way. At zero frequency, 
second-order poles contribute in a way 
which forces a careful consideration of the 
notation of summation and integration to 
avoid an error. A compact formula for the 
contribution of the second-order poles is 
presented. The sense in which the /-sum 
rule is satisfied is discussed. 
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1. The/Sum Rule and Transverse 
Response 

Causality — the principle that the state of a system 
depends upon its past but not its future — is expressed in 
the frequency domain as the principle that a causal 
response function has no poles in the closed upper half 
plane. Cauchy's theorem may then be used to derive the 
relationship between the real and imaginary parts of 
such a function, known in physics as the Kramers- 
Kronig relation. The dielectric function e{(0) is an 
example of such a function. 

Since e{(o) has no poles in the closed upper 
half plane, the same is true of [£(a>)- \]co. Cauchy's 
theorem implies that 

= hnjdo)[e(o))-l]o) (1) 



along the standard contour given in Fig. 1 . For a system 
of electrons in an electromagnetic field, the high- 
frequency response is given by the free electron form 



£ (0)) = l-% + O(a)- 4 ) 
or 



(2) 



where (O p is the plasma frequency given by (O^ = Ann^im, 
where n is the number density of electrons, m is the 
electron mass, and -e is the charge on the electron [1-3]. 
The OQ notation describes the limiting behavior [4]. 
We may break the contour integral into two parts and 
use Eq. (2) to simplify the integration along the semi- 
circle. Defining G) = Re l6 > 



= lim d(olme((o)(0 

£->« J-R 

2 

-f "d6Re ie -r^/te'* +0(R' 2 ). 
Jo /jV* 



(3) 
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Fig. 1. The contour of integration used to derive the /-sum rule. The limit of large radius, i.e., R^><*> will be 
taken. The contour of integration may be taken to be exactly on the real axis because technically, for any finite 
?7, the poles are in the lower half plane. The first-order poles at finite frequencies are shown with a small x. These 
arise from the fourth term in Eq. (9). The second-order pole at zero frequency is shown with a large x. 
Contributions to the second order pole arise from the second and third terms in Eq. (9). Note that these have the 
same small imaginary part, namely rj, which is related to the rate constant in the exponential which controls how 
fast the electric fields are turned on from a vanishingly small value in the distant past. 



The result 



K(o 2 p = J d(o\me{(o)(o 



(4) 



is known variously as the /-sum rule, the oscillator- 
strength sum rule, and the Thomas-Reiche-Kuhn sum 
rule. The /-sum rule may also be regarded as a high- 
frequency limit of the Kramers-Kronig relations [5]; 
in this case, the multiplication by ft) before Eq. (1) is 
better motivated. Physically the /sum rule puts an 
important constraint on the absorption of electro- 
magnetic radiation by a physical system. In this paper, 
we will consider how the /sum rule applies to the 
transverse response of the electron gas. 

The/sum rule should apply to all response functions 
that satisfy Eq. (2) in the high frequency limit. 
Although the exact dielectric function for an electronic 
gas is unknown, Lindhard [6] provides an analytic 
approximate form that incorporates the basic physics of 
the problem. Lindhard's dielectric function for the 
response of an electron gas to a longitudinal perturba- 
tion is widely discussed in textbooks [1-3], In his paper, 



Lindhard also presents the dielectric function for the 
response of an electron gas to a transverse perturbation 
[6], Little use has been made of the transverse response 
function, although, at least two articles report attempts 
to extend it to account for finite electron lifetime [7,8], 
and we used the function to discuss transverse response 
in the case of a fast electron traveling through the elec- 
tron gas [9]. Recently, current-current susceptibility for 
the electron gas has been derived as part of the study of 
dynamic exchange-correlation potentials [10,11]. The 
current-current susceptibility %jj is very closely relat- 
ed to the dielectric function (including both its trans- 
verse and longitudinal parts); specifically Eq. (4.165) 
of Pines and Nozieres [1] is (with a slight change of 
notation) 



CD 4 ne 2 
e(q,a>) = l — f- —XjA&oS) 

CD CD 



(5) 



where £ is the dielectric tensor, CD is the frequency, and 
q is the wave vector. In a basis which diagonalizes 
£ (q, CD), the dielectric tensor has two eigenvalues given 
by £%,©) and £ (1 \q,co) which are the transverse and 
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longitudinal dielectric functions, respectively [6,12]. 
The convention q = qq is used, where q is the magni- 
tude of the wave vector and q is a unit vector. 
Lindhard's transverse function may be obtained from 
the results of Bohm, Conti, Nifosi, and Tosi [10,11] 
by the application of Eq. (5). 

Lindhard gives the imaginary part of the transverse 
dielectric function as 



OK 



Im£ l ° (q, )) = ^ r x 



if}' 

— u(l-u 2 -z 2 ) for u + z<\ 
4 

— -[\-(u-z) 2 ] 2 for\u-z\<\<u + z 
12z 



(6) 







for I u - z I > 1 



for (O > 0. (An overall sign error has been corrected; 
odd parity applies for negative (O.) Lindhard defines the 
dimensionless variables z = qi(2k F ) and u = (D/(k F q), 
where (D = mcolh, h is Planck's constant, and k F is the 
Fermi wave vector. He gives the real part of the trans- 
verse dielectric function as 



Re£ (t \q,co)=l-^i-(z 2 +3u 2 +1) 

X[l-(z± W ) 2 ] 2 ln 



3 

32z 



z±u+\ 



z+u-\ 



(7) 



The sum over ± is just the two-term sum of the 
formula with ± —> + then ± —> -. This function 
obeys Eq. (2) independently for each q. It is even 
stated in the literature that the transverse dielectric 
function £ (0 satisfies the /-sum rule of Eq. (4) [13]. 

Let's check. The integrals are numerous but elemen- 
tary, and the result is 



\da)Jme (t) (q 9 (o)(0- 



710) 



^_l z 2 + _^ (1 _ z 2 )2ln 



3_ 
\6z 



\+z 



\-z 



(8) 



A single functional form exists for all z (physically 
all q) despite the fact that Eq. (6) has various analytic 
forms. Equation (8) is KG) 2 only for the limit of q -^ 0, 
whereas equality is expected for all q according to 
Martin [13]. The discrepancy is due to the existence of 
a pole term at Q)=0 [1,9]. Below, we rederive 
Lindhard's transverse dielectric function with an 
emphasis on the second-order pole. 



2. Transverse Response Revisited 

We specialize the derivation of Adler [12] for the 
dielectric function of a periodic solid with band theory 
to the case of the electron gas. We find 






2a) (k + q/2)®(k + q/2) 



nQ 2 



2^2 

m 0) p 



Ittt 



q-k + q 2 12 
{k + ql2)®{k + ql2) 



(9) 

■fv 



h z n ti(q-k+q 2 l2) 2 {±Q i + q-k + q 2 !2) 



where £1 is the frequency of the external perturbation 
including a small, positive imaginary part irj used in the 
"adiabatic turn-on" argument, i.e., £1 = (D+ irj, n is the 
number density of the electron gas, and scalars are 
promoted with identity tensors as required. Here, e is 
the dielectric tensor associated with the electric field 
derived from the vector potential, or equivalently, the 
dielectric function expressed in the Hamiltonian gauge 
which is defined by (/) = 0. The scaled frequency is 
defined by £l = m£l/ h. The /- are the Fermi occupa- 
tion numbers defined with Adler 's normalization con- 
vention, i.e., /- -21 F for a fully occupied state, includ- 
ing a factor of 2 for spin degeneracy; here V is the sys- 
tem volume. Adler's normalization convention leads to 



? / ^(2^) 3 



\dk, 



(10) 



as discussed, for example, in Ashcroft and Mermin 
[14]. For the electron gas, the energy of a state with 
momentum k is given by Z?- = h 2 k 2 /(2m), so 



E r - 

k+q 



■«-£M4- 



(ii) 



which appears in the derivation of Eq. (9). In deriving 
Eq. (9), we found that Eq. (5) needed to be corrected to 



£(q,co) = \- 



GK, 



Ane 2 



(co +irj) (co+iriy 



■XjA%co) (12) 



by the addition of irj to the two denominators. 



301 



Volume 113, Number 5, September-October 2008 

Journal of Research of the National Institute of Standards and Technology 



Equation (9) includes both longitudinal and trans- 
verse response. These are uncoupled [12]. The longi- 
tudinal response may be obtained by forming 
q-e(q,co)-q. In this case the second and third terms of 
Eq. (9) cancel, and the first and final terms give rise 
to Lindhard's longitudinal dielectric function, as we 
have verified in detail. (If the reader wishes to do so, 
we recommend forming the dimensionless variables 
K-qik F and w ± = 1/2 ± cd/q 2 to minimize the algebra). 
Because the longitudinal response is usually derived 
from the scalar potential, this equivalence is a conse- 
quence of gauge invariance which Adler has proved in 
the more general context of periodic potentials of solids 
within band theory [12]. 

To find the transverse response, the form <7± £(#>&>) 
■ q ± is required, where q ± is a unit vector orthogonal to q . 
The real part of Lindhard's transverse dielectric func- 
tion is given by evaluating Eq. (9). (The integrals are 
reasonably elementary and very similar to the longitu- 
dinal case; we made use of a partial fractions expansion 

v~ 2 (±€d+v)~ 1 = ±v~ 2 (d~ 1 -v~ l cd~ 2 +co~ 2 (±CD+vy l with 

v = qkx + q 2 12 and x = cos as well as \dtt n \n t = 



(/! + !> 



-Y +1 ln* 



(«+l) 2 f +1 . The imaginary part of 



Lindhard's transverse dielectric function reproduced in 
Eq. (6) is given by just the fourth term of Eq. (9). 

However, there remains an imaginary part associated 
with Qr 2 , which is defined to be (co+iJ])' 2 . In the 
Appendix, we derive the relation 



1 



(co + irj) 



= v ind{ co) _() 



(13) 



CD 



dco 



by an argument which is analogous to the argument 
used to derive the well-known relation [2] 



1 = P—-in8(aj). 

(co+irj) co 



(14) 



Here, V causes an integral to be evaluated as a princi- 
pal value, and 5 is the Dirac 5 function. 
Consider the dielectric function 



(01 



e(co) = l — p -, 
Q 2 



(15) 



Taking the imaginary part from Eq. (13), it is possi- 
ble to form the /-sum rule integral of Eq. (4) as: 



dCO(-CD 2 )(-x)8(CO) — CO = KCO 2 . ( 16 ) 

p dco p 



This example turns out to be the only relation that is 
not widely discussed which is needed to complete the 
discussion of the f-sum rule for Lindhard's transverse 
dielectric function. Equation (15) is the first and second 
terms of Eq. (9), so Eq. (16) represents their contribu- 
tions to the f-sum rule. (Of course, the first term, 1, 
has no imaginary part so has a zero contribution.) 
Equation (16) applies to the third terms of Eq. (9) as 
well, because the sum over k has no £1 dependence. 
Explicitly, the third term yields 



m I-- + -Z 2 -— (l-z 2 ) 2 ln 
Pi 8 8 16z 



1+z 



(17) 



which is just minus Eq. (8): the second term alone 
yields the expected answer of kco 2 . Alternatively, the 
pole term may be found by grouping the second and 
third terms. Application of Eq. (13) to Eq. (9) yield the 
pole term [9] 



W|2 + 2 z 2_J_ (1 _ z 2 ) 2l n 

lo 8 16z 



1+z 



(18) 



If this term is added to Eq. (8), the result is kco 2 . 

We emphasize that Eq. (8) and Eq. (18) were inde- 
pendently derived from Eq. (9) using Eq. (10), Eq. (13), 
and Eq. (14). We did not assume the value for the sum 
rule, but rather derived it. The discrete sum of Eq. (9) 
contributes a term at a pair of positive and negative 
frequencies, which are in general different for the 
different k and a zero frequency term for all k. As the 
number of k points becomes larger, the strengths in a 
small frequency interval differ by an ever-increasing 
factor, which accounts for the singular contribution at 
zero frequency. In other words, the contribution at zero 
frequency is a finite multiple of the contribution from a 
finite interval at non-zero frequency. Such a situation 
cannot be accounted for within elementary integral 
calculus. 

The traditional point of view is that the f-sum rule 
fails and the formula must be rewritten to include a pole 
term explicitly. Our point of view is that Eq. (18), when 
multiplied by the operator 5(co)d()/dco, should be con- 
sidered to be part of the imaginary part of £ (x \q,co), and 
that the f-sum rule is satisfied as written in Eq. (4). 
Because the imaginary parts of £ (t) (#, CO) arise from the 
same original formula, Eq. (9), and because the first- 
order poles and the second-order poles give analogous 
contributions, as shown in the Appendix, we see this ac- 
counting as more systematic, although less traditional. 
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3. Conclusions 

The key requirements for the^sum rule are causali- 
ty which implies there are no poles in the upper-half 
complex frequency plane, and the free-electron 
response at high frequencies. In the ^sum integrand, 
Im £(o))o), the real axis gets contributions from poles 
which are located just below it. First-order poles in 
£(a>) contribute for co^O, but a second-order pole in 
£(a>) contributes for a>=0. The formula for contribu- 
tions from such second-order poles is given by Eq. (13) 
which is analogous to the more familiar Eq. (14). It is 
somewhat unsettling that the principle "zero times any- 
thing is zero" does not apply to generalized functions as 
singular as those of Eq. (13). The traditional point of 
view is that the /-sum rule fails and a pole term is 
required. Our point of view is that the f-sum rule holds 
with a modified imaginary part of the dielectric func- 
tion. In physics, the Dirac 8-function is widely used, 
but it is just one of a number of generalized functions 
which have been studied for over 50 years [15]. By 
introducing a higher generalized function to the imagi- 
nary part of the response function, causality and the 
/■sum rule become more tightly linked. 

In the important special case of 



e(CD) = l- 



0)Z 



((D + ir\) 



(19) 



only the second-order pole contributes to the /-sum 
integral. The imaginary part of the dielectric function in 
Eq. (19) is rarely discussed because it is notationally 
awkward, without using notation such as that of 
Eq. (13), even though the/sum rule for the inverse of 
Eq. (19) is widely discussed as a plasmon. 

We appeal to authors and teachers to carefully distin- 
guish between co and Q = co+irj. The almost trivial 
difference between Eq. (5) and Eq. (12) took us sever- 
al weeks to fully comprehend. One of the big ideas of 
physics is that if the assumptions and rules of inference 
are stated clearly enough, verifying a derivation is 
purely mechanical. This only works if the notation is 
crystal clear. 



4. Appendix: Integral Relations for 

(co-\- irfy 1 and (a>+ irj)~ 2 

Equation (14) is widely used in theoretical physics 
[2]. The imaginary part of this relation may be demon- 
strated as follows for any smooth bounded function 
f((0): 



^ = limlm f dco^— 
^o J — (D + ii 



17] 



= -hm dco 7 \ rj 

77->oJ-» co +r\ 

defining A. Let u = co/rj. 

A= _ lim r du jm 

= -/(0)f" d »T^T 
= -nf(0) 



(A.1) 
(A.2) 



(A.3) 

(A.4) 
(A.5) 



establishing the result. The analogous derivation for 
(o)+/rj)" 2 is 

B = limlm f dco f ^ . (A.6) 

77->o J — ((Q + ir() 

= -21im f d<Q-4^%rm (A.7) 



„->oJ-~ (1+M 2 ) 2 ?] 




(A.8) 


- litnTiw " 2 f<n")-f(- 


-w) 


(A.9) 


— inn i ww _ _ 




do) l< "- OJ - (1+w 2 ) 2 




(A. 10) 


= -*/(0). 




(All) 
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Because the proof for the real part is nearly identical to 
helpful the well-established (g)+ /rj)" 1 case, it is omitted. This 

proves Eq. (13). 
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